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In order to explain how impurities affect the unconventional superconductivity, we study 
non-magnetic impurity effect on the transition temperature using on-site U Hubbard model 
within a fluctuation exchange (FLEX) approximation. We find that in appearance, the reduc- 
tion of Tc roughly coincides with the well-known Abrikosov-Gor'kov formula. This coincidence 
results from the cancellation between two effects; one is the reduction of attractive force due 
to randomness, and another is the reduction of the damping rate of quasi-particle arising from 
electron interaction. As another problem, we also study impurity effect on underdoped cuprate 
as the system showing pseudogap phenomena. To the aim, we adopt the pairing scenario for 
the pseudogap and discuss how pseudogap phenomena affect the reduction of Tc by impurities. 
We find that 'pseudogap breaking' by impurities plays the essential role in underdoped cuprate 
and suppresses the Tc reduction due to the superconducting (SC) fluctuation. 

KEYWORDS: high-Tc cuprates; impurity effect; fluctuation exchange; superconducting fluctuation; pseu- 
dogap 
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1. Introduction 

Until now many experiments on the effect of impurity 
substitution in hole-doped cuprate high-Tc superconduc- 
tors (HTS cuprates) have been done, and it is well known 
that the substitution in the copper site leads to a rapid Tc 
reduction. In their experiments, the following Abrikosov- 
Gor'kov (AG) formula"'^ is usually used to analyse the 
dependence of Tc on impurity content : 



^ ' 2 27rT, 



(1) 



where '4'{x) is the digamma function, Tco — T'c(nimp — 0), 
and r is the pair breaking scattering rate. In most of ex- 
periments which study impurity effect in the HTS, fitting 
the measured data by Eq.(l) seems to be successful, espe- 
cially in the low impurity concentration regime because 
of its linear bahavior. For example, it is reported theo- 
retically and experimentally that the Tc reduction by Zn 
substitution in HTS follows quantitatively Eq.(l), where 
Zn is treated as a strong non-magnetic scatterer, and F 
is evaluated in the unitarity limit. 

However, it is not so obvious that AG formula can be 
used to estimate impurity effects in HTS cuprates. The 
simplicity of Eq.(l) is based on the treatment of HTS 
as the weak-coupling d-wave BCS superconductor and 
the calculation technique of taking an average on Fermi 
surface. In such a treatment, the microscopic mechanism 
of superconducting pairing is not included at all, so the 
pairing interaction is not affected by impurities, whereas 
HTS is the typical system in which the pairing interac- 
tion is caused by Coulomb repulsion U and the momen- 
tum dependence plays the essential role in the mechanism 
of superconductivity. Thus, the effect of impurity on the 
pairing interaction may not be negligible in HTS. There- 
fore, it is necessary to study how impurities suppresses 
Tc from the microscpic view. To clarify the reason why 
Eq.(l) represents Tc reduction in HTS cuprate is one of 



the purpose of this paper and discussed in §3. 

As another problem, it is frequently debated that the 
impurity effect on underdoped cuprate is different from 
that on overdoped cuprate. In the experiments on the 
Tc-reduction by copper site substitution,^^' it has been 
noted that in underdoped regime there is more rapid re- 
duction of Tc than in the overdoped regime. In the other 
kind of experiments on electron irradiated samples, it is 
reported that the variation of Tc with defect content in 
underdoped and optimally doped YBa2Cu307_5 is quite 
linear and contrasts with the AG formula. Because one of 
the feature of underdoped cuprate in the clean samples is 
the pseudogap phenomena, it is natural to consider that 
such remarkable behaviors in underdoped region orig- 
inate from the pseudogap. Therefore, we disscuss how 
pseudogap phenomena affect the reduction of Tc by im- 
purities. This is the second purpose of this paper and is 
discussed in §4. 

Through this study, we treat impurities as non- 
magnetic ones in the mind of the Zn-substitution and 
defects; We assume that the scattering is elastic and that 
the impurity potential is the form of (5-function type and 
we take an average on the positions of impurities. More- 
over, we take the unitarity limit. 

This paper is constructed in the following way. In §2, 
we consider the impurity effect in the weak-coupling d- 
wave BCS model. We point out that Tc reduction in the 
model almost coincides with AG formula and that the 
slight deviation from the AG curve originates from the 
direct momentum summation done numerically. In §3, 
using the FLEX approximation, we discuss Tc reduction 
by impurities in the system where the pairing interaction 
is caused by U . First, we introduce the FLEX theory and 
our theorical treatment for impurities. In §3.2, we show 
how impurities affect the main factors to determine the 
value of Tc. In §3.3, we discuss the coincidence between 
the obtained Tc reduction and AG formula and we ex- 
plain the reason. In §3.4, we point out the necessity of 
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our treatment for the impurity effect in the HTS. In §4, 
using 'pairing scenario', we discuss how pseudogap phe- 
nomena affect Tc reduction due to impurities in the mind 
of the disordered underdoped cuprate. We first review 
the theorical treatment in the pairing scenario. In §4.2, 
we include impurity effect in the theory and show some 
features of Tc reduction in the presence of pseudogap. 
We point out that such features are caused by 'pseudo- 
gap breaking'. In §4.3, we explain what the 'pseudogap 
breaking' is. Finally, we shortly comment on the quan- 
tum dynamics of SC fluctuation. In the last section §5, 
we conclude the whole discussion. 

2. Reduction of Tc by impurities in the weak- 
coupling d-wave BCS model 

In this section, wc consider non-magnetic impurity ef- 
fect on Tc reduction in the unitarity limit within the 
weak-coupling rf-wave BCS model. We assume the fol- 
lowing separable attractive interaction: 

Vk,k' = -g^k^k', (2) 

where g > is the strength of the interaction and 
fk = coskx — COS ky, the form factor of (ix2_y2-wave. 
The transition temperature is determined by 

^ = -^J2^k,kG(k,tUJn)Gi-k,-ZUJn), (3) 
k,n 

where 

G(k, iuin) = [icOn - (eic - m) - Simp(ia;„)]~^ (4) 

with 

Simp {iuJn) = - 1 ^ "'"^ . T • (5) 

We numerically solve these equations. Here fj, is the chem- 
ical potential, a;„ is the fermionic Matubara frequency, 
Un = (2n + 1)ttT, and ni^p is the impurity concentra- 
tion. We consider the 2D square lattice and assume the 
following tight-binding dispersion 

Ck = — 2f(cos kx + cos ky) + At' cos kx cos ky. (6) 

The filling of electrons is represented by n, and n = 1 
corresponds to the half filling. The parameters are cho- 
sen as 2t = 1, t'/t = 0.25, g = 0.7t and n = 0.75 ~ 
0.90. The Fermi surfaces of bare electrons with disper- 
sion Eq. (6) and the density of state (DOS), p{w) = 
— i '^^lmG{k,Lo), are shown in Fig. 1. Here, w = is 
the Fermi level. 

When we assume the particle-hole symmetry and re- 
place the momentum summation with the Fermi sur- 
face average in calculating Eqs.(3)-(5), we obtain the 
well-known AG formula, Eq.(l), with the pair-breaking 
parameter F = -^j^-^'^ The critical impurity concen- 
tration n?^p is given by nf^^ = ^^-^0^, where 7' = 
1.781072. In the limiting case of low impurity concetra- 
tion, we have from Eq.(l), 

^ - 1 _ ^ r 

T ~ A ^ 1^ '^imp- [I) 

-tcO 4 riimp-lcO 

The numerically calculated results of Tc reduction and 
the fitted curves by Eq.(l) are shown in Fig. 2. The fit- 




' («0) ffl 

Fig. 1. Fermi surfaces of bare electrons with ek and DOS. The 

parameters are chosen as t' = 0.25t, g = 0.7t, T = O.OlSt and 
n = 0.75, 0.80, 0.85 and 0.90. 



ted AG curves are obtained in the following way. In low 
impurity concetration regime, the calculated reduction 
of Tc as the function of riimp is almost linear. Thus, we 
can estimate the coefficient of Eq.(7), ? — ^-^t^, which is 
independent of nimp, obtain the value of F and plot AG 
curves in each chart of Fig. 2. From the figure, we can 
see that Tc = Tc(nimp) obtained as the solution of Eq.(3) 
almost coincides with the AG curve. However, there is 
slight deviation of nf^^ from the value expected by AG 
formula. In the calculation of Eqs.(3)-(5), we have done 
directly the momentum summation not taking the aver- 
age on the Fermi surface. It results in the difference of 
the estimation of DOS and the difference leads to the de- 
viation of nf^p. In fact, from Fig. 1 and 2, we can see that 
as DOS increases at Fermi level, the deviation becomes 
large. Therefore, the deviation is the result of dispersion 
and filling of electrons. 

3. Reduction of Tc by impurities in the FLEX 
theory 

We start from on-site U t-t' Hubbard Hamiltonian as 
the model of the HTS cuprates. In order to describe the 
superconductivity in the Hubbard model with the repul- 
sive Coulomb interaction, it is necessary to derive the 
effective interaction from the many body effect. Here we 
adopt the FLEX approximation.^"^ Earlier calculations 
using the FLEX approximation were rather successful 
in determining many physical quantities in hole-doped 
cuprates. Therefore we consider the reduction of Tc by 
impurities within the FLEX theory. 

3.1 Model 

Here wc summarize the FLEX theory and explain our 
theorical treatment for impurities. 

We consider the t-t' Hubbard Hamiltonian on 2D 
square lattice. It is given by 

H = X](ek-M)cLcka+^ J2 4+qT4'-qiCk'iCkT, (8) 

kcr kk'q 

where cj^.^ is the creation operator of an electron with 
momentum k and spin a, ji is chemical potential, and U 
is the on-site Coulomb repulsion, ck is the same tight- 
binding dispersion as Eq.(6). The parameters are chosen 
as2t = l,t'/t = 0.25 and U/t = 3. These parameters are 
reasonable in the sense that the calculated Fermi sufaces 
coincide with those determined by the angle resolved 
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Fig. 2. Tc/Tco vs. nimp and the curves fitted by AG formula, Eq.(l). The fitting curves are obtained in the way explained in the text. 



n=0.90 
n=0.85 
n=0.80 
n=0.75 




Fig. 3. Fermi surfaces calculated by the FLEX approximation 
without impurities at T=0.015t. The change from the bare case 
is very slight. It is recognized by comparing this figure with Fig. 
1. 



photoemissioii (ARPES) experiments. Fig. 3 shows the 
Fermi surfaces calculated by the FLEX approximation 
without impurity at T=0.015t. The spectrum at (tt, 0) is 
above /i for n < 0.79. Comparing Fig. 3 with the ARPES 
data/" we can see that the above set of parameters cov- 
ers the optimal doped to overdoped region of the HTS 
cuprates, for example La2-xSrxCu04. 

In singlet pairing superconducting states, the one- 
particle Green's function is 2x2 matrix in Nambu no- 
taion, and expressed in terms of a self-energy matrix: 



G{k) = [G^\k)-t{k)]- 



where 



and 



iuJnTo - (Ck - IJ')t3, 



(9) 



(10) 



±ik) = iw„(l - Z{k))fo + X{k)T3 + cj){k)Ti. (11) 

Here A; is a shorthand notation as fc = (k, ia;„). Z{k) is 



the renormalization parameter, X(^k) is the energy shift, 
(p{k) is the gap parameter, and fj(i = 1 3) and fo 
are the Pauli matrices and the unit matrix, respectively. 
Z{k) and X{k) are defined as follows: 

Z{k) = l--^{i:{k)-E{-k)), (12) 



1 



X(fc) = -(E(fc)+I](-fc)), 



(13) 



where E(fc) is the normal self-energy. In these notations, 
the diagonal and off-diagonal one-particle Green's func- 
tions can be written as 

G{k) = -D{k){iuJnZ{k) + (ek - m) + X{k)), (14) 



F{k) = -D{k)m. 



with 



D{k) = [(a>„Z(fc))2 + ((ek-Ai) + X(fc))2 + 0(fc)2 



(15) 



(16) 



In the self-energy matrices, we simply neglect the 
mixed diagrams of Coulomb and impurity interaction: 



S(k,iw„) = SFLEx(k,ia;„) -|- Simp(ia;„). 



(17) 



The self-energy matrix Simp is given by Simp(«w„) = 
■fiimpH^Wn) where nimp is impurity concentration and t 
is t-matrix. The t-matrix is the solution of the following 
self-consistent equation: 



i{iLLln) =U + ^ G(k, iLJn)i{iuJn), 



(18) 



where u — ut^ and u is the amplitude of impurity poten- 
tial. In the unitarity limit, we obtain the diagonal and 
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the off-diagonal components of Simp as follows: 

5o(-«w„) 



Simp(iw„) = Hi 



™P5o(iw„)5(o(-iw„) +5i(ia;„)2' 



where 



and 



5o(«w„) = --^^G(k,ia;„), 

k 



(19) 

(20) 

(21) 

(22) 



SFLEx(k, iWn) is determined by the following spin and 
charge fluctuation mediated interactions: 



3 Xo(g) _ 1 s/ n' 



(23) 



(24) 



2 1 + UxUq) 2' 
Here q is a shorthand notation as g = (q, it/m), where 
= 2nTin is a bosonic Matsubara frequency. Xo and Xo 
are the irreducible spin and charge susceptibilities which 
are written as follows: 

xHq) = -^Y.iG{k + q)G{k) + F{k + q)F{k)], (25) 



kn 



kn 



In terms of these interactions, the diagonal and off'- 
diagonal components of the self-energy matrix, Sflex, 
are written as follows: 



SFLEx(fc) = ^ E Ves{q)G{k + q), 



qm 



with 



iFLEx(fc) = ^Y1 ^-(^ - ^')^(fc') 
k'.n' 



V,s{q)^Vs{q) + V,{q), 
VM = Vs{q)-VAq). 



(27) 
(28) 

(29) 
(30) 



Combining Eqs.(19), (20), (27) and (28), we obtain the 
explicit expression of Eq.(17): 

^(^) = mYl Ves{q)G{k + q) + Simp(ia;„), (31) 

qm 

m = Yl ^-(^ - k')D{k')ct>{k') + cf>,^p{iWn). 

k',7i' 

(32) 

Eq.(32) is the so called gap equation. We numerically 
solve Eqs. (9)-(32) in the self-consistent way, choosing the 



chemical potential so as to keep the flUing constant, 
" = 4 E '^•^ = E ^^n)e*""°+. (33) 



kn 



The divergent point of superconduncting susceptibil- 
ity, that is, the transition temperature T^, is estimated 
by the following eigenvalue equation: 

= E ^-(^ - k')D{k')U=^~^{k'). (34) 

^ k',n' 

The transition temperature is determined by the crite- 
rion that the maximum eigenvalue of Eq.(34) is unity 
at T = Tc. Here, we have assumed that </)(k, ia;„) is 
anisotropic, thus (i)imp(icj„) = 0, because in the clean 
case (j) has (ix2_y2 -symmetry. It is not immediately ob- 
vious that this assumption may be reasonable when the 
calculation is done self-consistcntly in the presence of im- 
purites. However, this is indeed the case as we will shortly 
disscuss in §3.5. Then, within the assumption, we study 
the impurity eff'ect on Tc obtained by solving Eq. (34). 

3.2 Impurity effects on the quasi-particle damping due 
to electron correlation and on the pairing interac- 
tion 

From the earlier calculation using the FLEX theory, 

it is well known that the main factors to determine 
the value of Tc are the magnitude of Ki(q, w) around 
q = (±7r, ±7r) and the imaginary part of self-energy 
around k = (±7r, 0), (0, ±7r). Therefore, we examine what 
is going on in these factors by impurities. 

In our calculation, Vs{q_, 0) has the commensurate peak 
at q = Q = (tt, tt) for n = 0.85 and the incommensu- 
rate peak for n = 0.75. As for the discussion below, the 
essential points do not change whether Vg has the com- 
mensurate peak or not. That is, the essential features 
discussed below do not depend on the amount of hole- 
doping. Then, we discuss mainly the commensurate case. 

Fig. 4 shows ImEpLEX at k = (0, tt) and ImT4(q, w) 
at q = Q as the function of w. The filling of electron is 
n = 0.85. In our formulation, the normal self-energy is 
given by S = SpLEX + '^imp- Then, llmSpLExl means 
the quasi-particle damping originated from the electron 
correlation. Prom the figure, we can see that and 
jlmSpLExl are suppressed by impurities. This is because 
the specific scattering process due to electron correlation 
is blurred by impurity scattering. By comparing the data 



of 'U = 2.71t, ni,„p = 0.00' and '17 : 



3i, niinp = 0.02', we 



can see that the amount of llmEpLEx] suppression by im- 
purities (riimp = 0.02) is almost identical to that by the 
reduction of U from 3t to 2.71t. In the usual FLEX the- 
ory without impurities, as U reduces, Tc becomes small. 
This is shown in Fig. 5. From the data of 'n = 0.85' 
in Fig. 5, the change of Tc expected by the suppression 
can be roughly estimated as ATc ~ 0.003t. It indicates 
that the reduction of Vs and llmSpLExl is rather large. 
Then, the pairing interaction and the pair-breaking due 
to the electron correlation are significantly suppressed by 
impurities. 
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Fig. 4. The imaginary part of Eplbx at k = (0, tt) (a) and 14 at Q = (tt, tt) (b) as the function of real lu. The filling of electron, n is 
0.85. The pure case corresponds to 'U = 3t, nij^p = 0.00' and the dirty case corresponds to 'U = 3t, nimp = 0.02'. For the comparison, 
the case of the reduced U without impuries {'U = 2.71t, nimp = 0.00') is included. The strength of the spin fluctuation interaction 
associated with momentum Q can be estimated by the integral, V^(Q, 0) = 2/^°° ^ImV^(Q, a))/a;. 




0.002 - 

I ' ' ' ' ' 

2.2 2.4 2.6 2.8 3 

U/t 

Fig. 5. Tc vs. U for various fillings of electrons without impurities. 
The feature that small U leads to small Tc is common in the usual 
FLEX theory. 

3.3 Reduction of Tc as the function of riimp and the 

mechanism 

Fig. 6 shows Tc vs. nimp obtained by the FLEX the- 
ory in the presence of impurities. We can see that the 
reduction of Tc ahuost coincides with AG curves and 
then keeps the almost same form as that in the weak- 
couphng case. This result is a surprise because in spite 
of the rather large pair-weakening and the reduction of 
the pair-breaking originating from the electron interac- 
tion by impurities as mensioned in the last section, those 
changes never seem to be effective on Tc. The result orig- 
inates from the cancellation between these two impurity 
effects. Fig. 7 shows Tc = Tc(niinp) in the FLEX theory 
with the impurity-suppressed pairing interaction, Vs, and 
the impurity-suppressed FLEX part of the self-energy. It 
can be seen that Tc hardly changes up to riimp ~ 0.03. 
This indicates that the cancellation between two effects 
occurs. 

Now, we summarize the Tc reduction by impurities in 
the FLEX theory. First of all, the value of Tc is mainly 
determined by the magnitude of 14 (Qj with q ~ Q 
and the quasi-pariticle damping rate |ImE(k, around 
k = (±7r, 0), (0, ±7r). K can be regarded as the factor 



to raise Tc and |ImS| as the factor to lower Tc. Impuri- 
ties significantly suppress both Vg and |ImS|; the pair- 
weakening and the reduction of the quasi-particlc damp- 
ing due to the electron correlation. These two effects, 
however, almost cancel each other to give no effect on 
Tc. Then, we can regard the system described by the 
FLEX theory with impurities as that with an impurity- 
independent pairing interaction and impurities; that is, 
the system can be described by the weak-coupling d-wave 
BCS model. Consequently, the reduction of Tc as the 
function of riimp keeps the almost same form as that in 
the BCS model. 

It is notable that the coincidence between the reduc- 
tion of Tc in the FLEX theory and AG formula is the 
very appearance. In order to analyze correctly the ex- 
perimental data on impurity effect in HTS, we must 
take into account the pair-weakening and the reduction 
of |ImScouiomb|, because the superconductivity emerges 
out of Coulomb interaction in HTS. In the next section, 
we will pick up the example which indicates that we have 
to take the reduction of |ImEcouiomb| into consideration. 

3.4 Dependence of impurity effect on k.-point 

We define 'hot spot' and 'cold spot' of the Fermi sur- 
face following Pines's definition for the hole-doping case. 
As we show in Fig. 8, 'hot spot' is located in the vicinity 
of the anti-ferromagnetic Brillouin zone, and 'cold spot' 
is far from the anti-ferromagnetic Brillouin zone. We find 
that there is a remarkable difference in the impurity ef- 
fect on the quasi-particle damping between hot spot and 
cold spot. At hot spot, the change of |ImI]| at Fermi level 
is small for the increase of impurity concentration, but at 
cold spot, |ImS| becomes large; that is, impurities hardly 
affect the quasi-particles at hot spot but are effective at 
cold spot. This is shown in Fig. 9. The feature of hot 
spot is Ek = Ck-Q = 0. Therefore, the quasi-particles at 
hot spot exchange the strongly enhanced spin fluctua- 
tion interaction Vg with the momentum Q, and have a 
large damping rate even in the absence of impurities. In 
the above discussion, we have found that impurities re- 





Fig. 9. The imaginary part of the quasi-particle self-energy in the presence of impurity scattering at hot spot (a) and cold spot (b). 
The parameters are chosen as n = 0.85, t' = 0.25t, U = M, and T = O.OlSt. 



duce the spin fluctuation part of |ImE|, i.e. llmSpLExl- 
Therefore, at hot spot, the spin fluctuation part and the 
impurity part of |IniE| preserve a balance, and there- 
fore, ImS = ImSpLEX + ImSimp hardly change. On the 
other hand, at cold spot, since the original quasi-paritcle 
damping rate by the spin fluctuation scattering in the ab- 
sence of impurities is small, the reduction of llmEpLEx] 
is very slight, and then impurity scattering makes the 
magnitude of |ImS| large. 

The following point should be noted. In the absence 
of impurities, the conductance is mainly determined by 
the contribution from cold spot region. As impurities in- 
creases, the momentum dependence of lifetime becomes 
small, and the conductance becomes to be determined 
by the whole k's on Fermi surface. 



3.5 The supplementary discussion 

Here, we make sure that the assumption stated at the 
last of §3.1. Fig. 10 shows the momentum dependence of 
the gap parameter for the clean case and the dirty case. 
In both cases, the dx^-y^-wave arises clearly. Therefore, 
within the model we explained in §3.1, impurities do not 
change the symmetry of the superconducting gap. Then, 
the neglect of 4>{i-a-p in Eq.(34) is justified. 

4. Impurity effect in the presence of the pseu- 
dogap 

Pseudogap phenomena can be seen in underdoped 
cupratc. In order to describe the pseudogap, we adopt 
the "pairing scenario" where the pseudogap is induced 
by the superconducting (SC) fluctuation. ^^^^^ In the ap- 
proach, pseudogap phenomena are understood as the res- 
onance phenomena between the states of Cooper pairs 
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Fig. 7. 

theory with impurity-suppressed Sflex and U = 3t without 
and with Simp, respectively. The filling of electrons is 0.85. Dot- 
ted line is as same as the data of 'n = 0.85' chart in Fig. 6, 
although here Tc is not scaled by Tco(= 0.0137t). This is shown 
for compajrision. 



(0,7[) 



cold spot- — /\ 






^\ hot spot 



Fig. 8. The definition of 'hot spot' and 'cold spot' of the Fermi 

surface. The Fermi surface is obtained for the parameters chosen 
as n = 0.85, t' = 0.25, !7 = 3t and T = 0.015t. 



and quasi-particles in the Fermi liquid; the resonance is 
described by extraordinarily increased imaginary part of 
the quasi-particle self-energy. Y. Yanase and one of the 
authors (K. Y) have succeeded in describing the pseu- 
dogap by the 1-loop order theory with respect to SC 
fluctuation on the basis of the Hubbard model with d- 
wave attractive interaction.^^ Therefore, including effects 
of impurities into the theory, we discuss the Tc suppres- 
sion by impurities in the presence of the pseudogap. It 
is notice that the similar study has already been done 
by Qijin Chen and J. R. Schrieffer.^" However, their ap- 
proach for pseudogap is different from ours. They have 
used a GqG approximation scheme. Thus our study may 
serve as a test to judge what is an appropriate theory of 
pseudogap. 



Fig. 10. The momentum dependence of the gap paremeter (/'(k, 0) 
in the superconducting state. The parameters are chosen as t = 
0.5, t' = 0.25t, T = 0.007t, U = it, and n = 0.85. 




(^0) 



Fig. 11. The Fermi surface calculated using the parameters cho- 
sen as t' = 0.25t and n = 0.9. The point A shows k = (tt, 0.127r) 
and B shows k = (0.567r, O.SItt). A is so called antinodal point. 
We refer these points A and B many times in the below discus- 
sion. 



4-1 Reviews on the 1-loop order theory and the self- 
consistent calculation 
The attractive Hubbard model on 2D square lattice is 
given by 



kcr k,k' q 



Cq_k|Cq_k'iCk'T- 

(35) 

Here, Vk.k' is the rfx^-y^-wave separable pairing interac- 
tion written by the same form as Eq.(2) and ^k = Cfc — 
where ek is the same dispersion as Eq.(6). The parame- 
ters are chosen as t' ~ 0.25t, g = 0.4. the filling number 
n = 0.9, and t = 0.1 ~ 1.0. Fig. 11 shows the Fermi sur- 
face. The physical meaning of changing the value of t is 
as follows. Generally, in the high density electron system, 
the band width W (= 8t) of quasi-particle is renormal- 
ized by the effect of the electron correlation. Then, small 
t corresponds to the high density case. As for hole-doped 
HTS cuprates, small t corresponds to the underdoped re- 
gion and large t does to the overdoped region. Therefore, 
by changing the value of t, we can cover the whole region 
of cuprate. 

In the 1-loop order theory, the self-energy by SC fluc- 
tuation is given by the diagram in Fig. 12 and written 



8 J. Phys. Soc. Jpn. 



Full Paper 



Kiwamu Kudo and Kosaku Yamada 




Fig. 12. Sclf-cncrgy diagrams in the 1-loop order theory. The dot- 
ted line means k' • 



as 



I]sc(k, iiOn) = XI *^»»)G(q - k, Wrn - ii^n), 

(36) 



with 



-9 



1 -flXpplq.^^^m) 

is the bosonic Matsubara frequency, 



(37) 



Here v„ 

2iTTm, and Xpp is the pairing susceptibiUty written as 

Xpp(q, i^m) = Jp^Yl "^kG^C*' i^n)G{<l - k, iVm - iuJn)- 

(38) 

In the 1-loop order theory, the Green functions, 
G(k, zcij„), in Eq. (36)^(38) are given by bare ones, 
Go(k,it^„) = {iu)ri -^k)"^- 

Fig. 13 shows the real and imaginary parts of the self- 
energy and the spectral weight obtained by the 1-loop 
order theory. It can be seen that as a result of the reso- 
nance scattering, the damping rate |ImE^(,(k, w)| at the 
Fermi level is large and then the spectral weight is sup- 
pressed there. The suppression of the spectral weight is 
related to the suppression of DOS. Hence, as mentioned 
above, the 1-loop order theory well describes the gap-like 
behavior of the single-particle spectrum, that is the pseu- 
dogap phenomena. It also can be seen that ReE^(k, w) 
has a positive slope at the Fermi level. It is notable that 
this non-Fermi liquid behavior is naturally derived from 
the Fermi liquid state. 

It should be paid attention that the above calucula- 
tion by the 1-loop order theory is applicable to the higher 
temperature than the BCS mean field transition temper- 
ature, . In the framework here, the transition tem- 
perature is determined by the so called Thouless condi- 
tion: 

l-,9Xpp(O,0) = 0. (39) 

For the bare Green function Go(k, iw^), Eq.(39) is iden- 
tical with the BCS mean field form: 

. 1 

' j3N 



k,n 



(40) 



Therefore the transition temperature in the 1-loop order 
theory is identical to T^'^ . As the result, the temperature 
region as to the pseudogap phenomena is T > . 

To include the suppresion of Tc by the effect of the 
SC fluctuation, we have to go beyond the 1-loop order 
theory . In order to treat such effect, we carry out the self- 
consistent calculation (the self-consistent 1-Ioop theory) 
using the renormalized Green function, i.e., = Gq ^ — 
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Fig. 14. t vs. Tc in the self-consistent 1-loop theory and the BCS 
mean field theory. 'BCS' means the solution of Eq.(40), that is 
jAiF rpj^g parameters are chosen as g = 0.4 and n = 0.9. 



In the self-consistent 1-loop theory, we solve 
Eq.(36)~(38) self-consistently. choosing the chemical po- 
tential n so as to keep the filling constant n = 0.9. Using 
the self-consistent solution, we estimate Tc by the Thou- 
less condition. It should be noticed that in 2D system 
the superconducting transition does not occur as is well 
known as Marmin- Wagner's theorem. Therefore, we phe- 
nomenologically introduce the three-dimensionality and 
define Tc as the temperature in which 1 — 5Xpp(0)0) = 
0.015. By the procedure, we obtain the phase diagram. 
Fig. 14, which shows t vs. Tc in the self-consistent 1-loop 
theory and the BCS mean field theory. In the figure, as 
t decreases, Tc is more; suppressed from the value of the 
BCS transition temperature. The small t leads to large 
SC fluctuation via large g/t\ therefore Fig. 14 shows the 
suppression of Tc by the effect of SC fluctuation. This 
suppression is mainly due to the reduced DOS. Actu- 
ally, the smaller is t, the more suppressed is the spectral 
weight at Fermi level. This is displayed in Fig. 15 which 
shows the real and imaginary parts of the self-energy and 
the spectral weight obtained by the self-consistent 1-loop 
theory. It can be seen that at t = 0.05, |ImE^(k, a;)| is 
most enhanced and then the spectral weight near the 
Fermi level is most suppressed. Unlike in the 1-loop or- 
der theory, the Fermi liquid behaviors recover as can be 
seen in the negative slope of ReS^(k, w) at the Fermi 
level. Moreover, the resonance feature is weaken and the 
spectral peak at Fermi level is somewhat recovered. Then 
the gap-like behavior disappears. This is due to the par- 
tial summation included in the self-consistent calcula- 
tion. However, the reduction of DOS at Fermi level as the 
feature of the pseudogap remaines. Yanase^^ has shown 
that the peak structure near the Fermi energy in the 
spectral weight is not reliable but Tc obtained by the self- 
consistent 1-loop theory is reliable, since the value of Tc 
is determined by wide structure in the spectrum. There- 
fore, we consider the impurity effect on Tc in the presence 
of the pseudogap within the self-consistent treatment. 
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Fig. 15. The real and imaginary parts of the self-energy and the spectral weight obtained by the self-consistent 1-loop theory. Here, 
g = 0.4 and T = 0.02. The hopping parameter t is chosen as t = 0.05, 0.1 and 0.2. 




Fig. 16. Impurity contributions to Xpp which are formally taJsen 
into a<;count in our calculation. 

4-2 Self- consistent 1-loop theory with impurities 

Let us include impurity effects in the self-consistent 
1-loop theory. We incorporate the impurity effect into 
the single-particle Green function via the self-energy and 
also take into account the diagrams shown in Fig. 16. 
However, in our d-wave model, the diagrams in Fig. 16 
have no contribution due to the momentum summation 
which has the form of -i- (/^^ • • • , where (ySk = cos — 
cos ky. Therefore, impurities enter just only the Green 
function via Dyson equation: 

G(k, iujn)-^ = GoiK iwn)~^ - S(k, iWn), (41) 

with 

S (k, iw„ ) = Ssc (k, iw„ ) + Simp (i w„ ) . (42) 

Here 

We solve Eqs.(36), (37), (38), (41), (42) and (43) self- 
consistently choosing /z so as to keep the filling of elec- 
trons. 

Fig. 17 shows nimp vs. in the self-consistent 1-loop 
theory and, for comparison, that in the BCS model. In 
the self-consistent 1-loop theory, the critical impurity 
concentration for t = 0.2 is larger than that for t = 0.3 




0.05 0.1 0.15 0.2 
Himp 



Fig. 17. nimp vs. Tc in the BCS model and in the self-consistent 
1-loop theory for the case of t = 0.2 and 0.3. The parameters are 
chosen as g = 0.4 and n = 0.9. 

in spite of smaller Tco for t = 0.2 than for t — 0.3. This 
situation is contrary to that in the BCS model without 
SC fluctuation. In §2, we have seen that Tc reduction 
in the BCS model almost coincides with AG curves and 
that nP^p is given by ^2°]^°° in AG formula, that is, 
n?jjp depends on T^o. Therefore, we can consider that 
the magnitude of n^^p in the self-consistent 1-loop the- 
ory is mostly determined by the mean-field value of n-j^p 
which reflects the value of . Strictly speaking, how- 
ever, there exists the clear difference between n\^^ in the 
mean-fleld theory and that in the self-consistent 1-loop 
theory, as seen in Fig. 17. We consider that this difference 
originates from a quantum dynamics and will comment 
on this issue in §4.4. 
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Fig. 18. The Tq scaled by Tco as the function of the impurity 
concentration. The parameters are chosen as g = 0.4, n = 0.9 
and t = 0.1, 0.2 and 0.3. 
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Fig. 19. The xiimp vs. Tq data in the self-consistent 1-loop theory 
for the parameter as t = 0.2, g = 0.4, and n = 0.9, and the 
fitting curve by Eq.(l). The fitting is done in the same way as in 
fi2. 



The most peculiar feature of reduction is its nearly 
linear variation with impurity concentration. This is 
shown in Fig. 18 and 19. From Fig. 18, we can see that 
the curvature of Tc = rc(ninip) gradually becomes posi- 
tive as t is small and the reduction of is almost linear 
at t = 0.2. This feature originates from the strong SC 
fluctuation. In Fig. 19, we compare the data for t = 0.2 
and AG curve. AG formula can not be the good fitting 
curve due to the absence of the negative curvature in 
the calculated data for t = 0.2. Rullier-Albenque et al. 
have reported this kind of the linear variation. ^"^ They 
have measured Tc and ab plane resistivity of electron 
irradiated YBa2Cu307_5 crystals (underdoped and op- 
timally doped) and mensioned that the variation of Tc 
does not follow any current prediction of pair-beraking 
theories. However, we believe that a pair-breaking theory 
still holds even in their experiment and the 'pseudogap 
breaking' by impurities or defects, which will be intro- 
duced below, is the reasonable keyword for an explana- 
tion of their data. 

So, why the linear variation of Tc seen in Fig. 19 
emerges? 

Before tackling the issue, we analyze carefully Fig. 20 
and Fig. 21 which show the spectral weight and the imag- 
inary part of S and T,sc at A and B in the k-space as the 



funtion of uj for nimp = 0.00, 0.07, 0.14. We proceed with 
a discussion regarding the point A as the representative 
of the k-points around the 'hot spot' and the point B as 
the one around the 'cold spot', and then simply call A 
'hot spot' and B 'cold spot' in the below discussion. 

Comparing Fig. 20(a) and Fig. 21(a) for nimp — 0.00, 
we find that the spectral weight at hot spot has a short 
peak but the one at cold spot has a sharp and high peak. 
This difference owes to the form factor which has a 
large value at hot spot and a small value at cold spot. 
Physically, it can be seen that because of SC fluctuation, 
there is almost no coherence as the well-defined quasi- 
particle at hot spot even in the clean case. Let us in- 
troduce impurities. We find the remarkable feature that 
whereas the spectral peak at cold spot is largely sup- 
pressed by impurities, the one at hot spot is hardly af- 
fected by impurities and keeps the almost same form as 
that in the clean case. Owing to the direct connection 
between the spectral weight and the imaginary part of 
the self-energy, this kind of feature also can be seen in 
Fig. 20(b) and Fig. 21(b) which show the quasi-particle 
damping, |IniS|. As impurity concentration increases, 
while the lifetime of the quasi-particle at cold spot be- 
comes short, the one at hot spot hardly changes. That 
is, impurities seem to mainly affect the quasi-particles 
at cold spot and have no contribution to the ones at 
hot spot, although the contribution of impurities to the 
quasi-particle damping should not depend on momen- 
tum k. The reason that such things arise is clarified by 
taking accout of Fig. 20(c) and Fig. 21(c) which shows 
ImEsc, the contribution of SC fluctuation to the quasi- 
particle damping. These figures have the common feature 
that jlmSscI becomes small by impurities. It means that 
as impurities increase, the quasi-particles tend not to be 
affected by the SC fluctuation. In terms of OTir formula- 
tion, this is due to the diagram such as Fig. 22, which 
arises from the process of the self-consistent calculation. 
Impurities interfere in the scattering by the SC fluctua- 
tion. Prom these understandings, we can now appreciate 
the behavior that |IniS[ at hot spot keeps the almost 
constant value at Fermi level. This is because the con- 
tribution of ImE^ 
the cancellation makes ImS 



cancels out the one of ImSimp, and 



ImS, 



ImEi, 



^inip nearly 

constant. At cold spot, due to the smallness of the sup- 
pression of |ImSsc|, the cancellation is not complete, and 
impurities effectively contribute there. 

The idea of the suppression of llmEgd allows us to 
tackle the above issue. On the basis of the idea, we have 
drawn Fig. 23 which gives us a rough schematic under- 
standing for the linear variation of Tc. From the figure, 
the value of Tc in the self-consistent 1-loop theory can 
be expressed by the following form: 

jiMF 



where STw 



■5Ti-{STi^p-ST2), (44) 

etc. are explained in the caption of Fig. 23. 



As impurity concentration increases, not only STi^p but 
also 6T2 increases via the suppression of |ImEsc|- Then, if 
<^Timp — (ST2 were in proportion to nimp, Tc would decrease 
linearly according as nimp picks up. In this sense, our 
parameters used to draw Fig. 19 is such ones that keep 
the delicate balance between 6T2 and JTimp which leads 
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Fig. 22. The example of the self-energy diagram included in the 
self-consistent calculation. The wavy line means SC fluctuation 

interaction and the dotted lines with the cross mark the impu- 
rities scattering averaged over on their positions. This type of 
diagram makes weak the contribution of SC fluctuation to the 
self-energy. 
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Fig. 23. The figure for the rough understanding that To is sup- 
pressed almost linearly. ST^^p means the magnitude of the Tc 
reduction form T^^ in the BCS mean-fleld theory. STi means 
the difference of T^^^ and Tco in the self-consistent 1-loop the- 
ory, that is, the magnitude of suppression by SC fluctuation 
without impurities. The 'adjoint line' is drawn by shifting the 
'BCS(t = 0.2)' data down by the amount of STi. ST2 means the 
decline in the the magnitude of Tc reduction by SC fluctuation. 



to the nearly linear variation of Tc as the function of 

fiiuip- This is the answer for the above issue. 

The essential and important point of the discussion up 
to here was the idea that |ImEsc| is suppressed by im- 
purities. We call the idea 'psuedogap breaking'. We may 
doubt that the word is not consistent with the idea be- 
cause there is no gap-like behavior in the spectral weight 
of Fig. 20 or Fig. 21. However, the pseudogap is indeed 
destructed by impurities as we will explain in §4.3. 

Another important result we have to point out is 
included in Fig. 24 which shows the comparison of 
the quasi-particle damping by impurities in the self- 
consistent 1-loop theory and that in the mean-field the- 
ory. From the figure, we can see that the former is larger 
than the latter. As we mentioned in §4.1, SC fluctuation 
results in the reduced DOS. In our formulation, impuri- 
ties arc treated in the unitarity limit, then the reduced 
DOS loads to a large ImSimp as we can understand by 
Eq.(43). Therefore, the behavior shown in Fig. 24 can 
be intcrploted by the reasoning that the contribution of 
impurities to the quasi-particle damping becomes effec- 
tively large by the effect of the SC fluctuation. 

4-3 Effective 1-loop model 

In this section, we explain the concept of the 
'pseudogap-breaking' by impurities. To the aim, we in- 



Fig. 24. The impurities contribution to the imaginary part of self- 
energy as the function of u in the self-consistent 1-loop the- 

ory(SlL) and in the mean-field theory(BCS). The parameters 
are chosen as t = 0.2, g = 0.4, n = 0.9, and T = 0.02. 



troduce an effective model for the self-consistent 1-loop 
theory with impurities. The feature of this model is the 
maintenance of the gap-like behavior as can be seen in 
the 1-loop order theory even though the calculation to 
determine Tc is done self-consistently. 

4.3.1 Derivation of the effective 1-loop model 

Since T{Q) has the sharp peak around Q = (q, ii^m) = 
for T ^ Tc, we can expand TiQy^ around Q = 0. 
Moreover, we use the static approximation, i.e. we only 
take into account i'm = ^ for T{ci,iy.m). Then 



T(q,0) 



1 - .9Xpp(q,o) 

can be written as the following form: 

-9 



T(q,0)^ 



r + hc0 ' 



where 6 is a constant and 

r-=l-5Xpp(0,0). 



(45) 



(46) 



(47) 



By making use of Eq.(46), we can rewrite Eq.(36) as 

follows: 



Ssc(k, iujn) 



-9 



I3N ^ r + 6q2 



-lUlr, 



Sq-k 



-'imp(~*'A'n) 



-,.2 9 



-y 



q 



1 



«Wn + + Vxqx + Vyqy -\- Siinp(-iWn) 

(48) 

Here, ^ = b/r, and we have replaced with 

t-l <r-l 

^ , 2 c-2 I 2 ■ Moreover, we have used that £a-k around 



q = 



can be written as 

$q-k = ?k 



(49) 



12 J. Phys. Soc. Jpn 
(a) 

10 



Full Paper 



Kiwamu Kudo and Kosaku Yamada 







nimp=0.00 

nimp=0.07 

nimp=0.14 




/ 










-0.4 


-0.2 


} 0.2 0.4 




Fig. 20. The spectral weight and the imaginary part of S and Sgc as the function of real frequency u) obtained by the self-consistent 
1-loop theory with impurities. Here, t = 0.2, g = 0.4, T = 0.02, and k is at A in Fig. 11. The impurity concentrations nimp are chosen 
as nimo = 0.00, 0.07 and 0.14. 
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Fig. 21. The spectral weight and the imaginary part of S and Sgc as the function of real frequency ui obtained by the self-consistent 
1-loop theory with impurities. Here, t = 0.2, g = 0.4, T = 0.02, and k is at B in Fig. 11. The impurity concentrations riij„p are chosen 
as riimp = 0.00, 0.07 and 0.14. 



where 



and 



Ua;(k) = —2smkx{t — 2t' cos ky), 



(50) 



Vy{k) = —2smky{t — 2t' coskx), (51) 
and defined as follows: 



A^ = -^En<,,o). 

q 

By using this A^, we can express ^ as 



4/3r A2 ■ 



(52) 



(53) 



Integrating the last expression of Eq.(48), we finally ob- 
tain 



ioJn + + i{\Vx\ + \vy\)i ^SgnWn -I- Siinp(-iw„ 



(54) 



(55) 



The renormalized Green function is given by 

G{]&,iujn) = -. T ^ ^. — T — T. ^ 

- ?k - Ssc(k,za;„) - S;imp(«w„) 

where Sinip(ia'„) is as same as Eq.(43). 

We solve Eqs.(38), (45), (52), (53), (54), (43) and (55) 
self-consistently choosing the chemical potential ji so as 
to keep n = 0.9. 



4-3.2 Results of calculation and 'pseudogap breaking' 

Fig. 25 shows t-T^ phase diagram in the effective 1-loop 
model without impurities and for comparison that in the 
self-consistent 1-loop theory we have already displayed 
in Fig. 14. In the effective 1-loop theory, the larger sup- 
pression of Tc from the BCS value than that in the self- 
consistent 1-loop theory can be seen. That may be be- 
cause wc neglect the components of 7^ in r(q, ifm) 
and then the SC fluctuation as a thermal fluctuation is 
large in the effective model. 

Fig. 26 shows nimp-Tc phase diagram in the effective 
1-loop theory with the parameter is the same as in Fig. 
17. The rough behavior of is identical to that in the 
self-consistent 1-loop theory. 

The obvious difference between the effective model and 
the self-consistent 1-loop theory is reflected in the spec- 
tral peak at hot spot which is shown in Fig. 27. The 
gap-like behavior in the spectral weight, that is, the pseu- 
dogap, can be seen in the data for m 



imp 



0.00. It is no- 
table that as riimp increases, jlmSsd becomes small, the 
pseudogap is suppressed, and then the single peak in the 
spectral weight emerges. This is nothing but 'pseudogap 
breaking'. Now, we have understood that the suppression 
of jlmEscI by impurities is closely related to the pseudo- 
gap breaking. 

As is already mensioned, the suppression of jlmSsd 
i.e. the pseudogap breaking, results from the self-energy 
diagram with such a form as Fig. 22. Because the causal 
relationship is very clear in the effective 1-loop theory, 
let us consider the relation again. 
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Fig. 25. t vs. Tc in the self-consistent 1-loop theory and in the 
effective 1-loop model. The parameters are chosen as g = 0.4 
and n = 0.9. The data of 'self-consistent 1-loop' is as same as 
the data in Fig. 14. 
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Fig. 26. 

^imp Tc in the effective 1-loop model for the case 
of t = 0.2 and 0.3. The parameters are chosen as g = 0.4 and 
n = 0.9. 



For simplicity, we adopt the following form of S 
Si„ip(«w„) = -i^sgnoJn- 



imp) 



(56) 



Here, t(oc ^j^^p) is the lifetime of the quasi-particle 
by impurity scttering. By using this expression, we can 
rewrite Eq.(54) as 



i(^n + + irkSgna;„ ' 



(57) 



where 



rk = 



1 

+ 27- 



(58) 



It is notable that the expression of Eq.(57) is nothing but 
the diagram of Fig. 22 itself within the approximation by 
which we have used to derive the effective 1-loop model. 
By analitical continuation, iuin co + i5, in Eq.(57), we 
obtain 

AV^rk 



Iml]sc(k,a;) = - 



(59) 



From this equation, we find that |ImSsc(k, a;)| has a peak 
at w = for the k on the Fermi surface and its top value 

^2 2 

is -Y^- Therefore, as riimp increases, Fk becomes large, 
and then the peak of |ImEsc(k, w) | lowers. This tendency 
coincides with the behavior of ImEgc in Fig. 27. As the 
result, we recognize that the scattering process such as 
Fig. 22 is essential for pseudogap breaking. That is, the 
broadening of the one-particle state by impurities makes 
weak the scattering due to SC fluctuation and it leads to 
the pseudogap breaking. 

4-4 Comments on the quantum dynamics of SC fluctu- 
ation 

In §4.2, we have mensioncd that the difference between 
n-jjjp in the BCS theory and that in the self- consistent 
1-loop theory. If the SC fluctuation were perfectly ther- 
mal, rtfujp's in the both theories would coincide with each 
other, because the thermal fluctuation should vanish at 
T = 0. In the self-consistent 1-loop theory, however, the 
quantum dynamics of SC fluctuation is included. Actu- 
ally, in Eq.(36), all i^m's are included in the summation 
not just only = 0. Therefore, we expect that the 
above difference originates from the quantum dynamics. 
The expectation is supported by Fig. 28 which shows 
"■imp vs. Tc in some theories. The 'staticSlL' in the fig- 
ure means that we have applied the static approximation 
to Eq.(36), that is, taken just only f „ = into account 
in the self-consistent 1-loop theory. In Fig. 28, it can 
be seen that n'^^^ in the three theories except the self- 
consistent 1-loop theory almost coincide with each other 
and n?jjp ~ 0.21. The reason would be that in 'staticSlL' 
and 'effective', only the thermal fluctuation is taken into 
account. On the other hand, nf^^ in the self-consistent 
1-loop theory is smaller than that in the other theories 
and nfjjjp ~ 0.15. Then, we consider that the quantum 
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Fig. 28. 

^imp Tc in tliG self-consistent 1-loop thcory(SlL), SIL 
witliin tlic static approximation(staticSlL), the cfEcctive 1-loop 
modcl(cfl'cctivc) and the BCS mean-field tlicory(BCS). The pa- 
rameters are chosen as t = 0.2, g = 0.4 and n = 0.9 in all 
theories. 



dynamics of SC fluctuation plays a role to reduce n'(^^. 
Farther understanding about the quantum fluctuation is 
one of our future problems. 

5. Conclusion 

The reduction of Tc due to impurities in the FLEX the- 
ory is almost same as that in the weak-coupling d-wave 
BCS model. This is because the cancellation between two 
effects by impurities; one is the pair-weakening, and an- 
other is the reduction of the pair-breaking produced by 
electron correlation. Therefore, even in the HTS, where 
the pairing interaction is caused by U , AG formula still 
can be used for rough estimation of reduction due to 
impurities. However, the following points should be no- 
ticed. By the numerical calculation, we have found the 
deviation of the critical impurity concentration from the 
nfjjp expected by AG formula. It depends on the flne 
structure of the dispersion and the Fermi surface via 
DOS. The larger is DOS, the larger the deviation of n\^^ 
is. In order to analyze the hole-doping dependence of Tc 
reduction by impurites in detail, we should notice the 
form of Fermi surface not just only adopting AG for- 



mula. 

As for the dependence of impurity effect on k-point, 
impurities mainly affect the quasi-particles at cold spot. 
This is because impurities reduce the inelastic scatter- 
ing by spin fluctuation which is effective to the quasi- 
particles at hot spot. This is the good example to indi- 
cate that we ought to take into account the reduction of 
the damping rate of quasi-particle arising from electron 
interaction. 

In §4, we have studied how pseudogap phenomena af- 
fect the reduction of due to impurities. The varia- 
tion of Tc as the function of nimp is too linear to fit by 
AG formula in the region of rather large SC fluctuation. 
The linear variation is caused by the suppression of the 
quasi-particle damping arising from SC fluctuation i.e. 
the pseudogap breaking. 
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